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Abstract 

We calculate the nonlinear magneto-optical response of free-standing fee 
(001), (110) and (111) oriented Fe monolayers. The bandstructures are de- 
termined from first principles using a full-potential LAPW method with the 
additional implementation of spin-orbit coupling. The variation of the spin- 
orbit coupling strength and the nonlinear magneto-optical spectra upon layer 
orientation are investigated. We find characteristic differences which indicate 
an enhanced sensitivity of nonlinear magneto-optics to surface orientation and 
variation of the in-plane lattice constants. In particular the crossover from 
onedimensional stripe structures to twodimensional films of (111) layers ex- 
hibits a clean signature in the nonlinear Kerr-spectra and demonstrates the 
versatility of nonlinear magneto-optics as a tool for in situ thin-film analysis. 

75.30.Pd;78.20.Ls;73.20.At;75.50.Bb 
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I. INTRODUCTION 



The Nonlinear Magneto-Optical Kerr effect (NOLIMOKE) is a unique optical tool to 
analyze tliin-film systems, since it is exclusively generated at surfaces and interfaces, where 
the local inversion symmetry is broken Therefore is has attracted considerable interest 
in recent research on interface magnetism |[^-[T0|. In this paper we investigate the depen- 



dence of the nonlinear magneto-optical Kerr-spectra and their microscopic origin [|Tl|] 
the spin-orbit coupling (SOC) and the magnetic moment - of Fe monolayers on structural 
changes. In particular the effects of different in-plane lattice constants, different coordina- 
tion and onedimensional stripe structures are studied. The results were obtained using the 
full potential linearized augmented plane wave (FLAPW) method WIEN95 [l^ with the 
additionally performed implementation of spin-orbit coupling. 

Two features are responsible for the strong interest in magnetic thin films: (i) the mag- 
netic properties strongly depend on structural changes and (ii) the spin-orbit induced effects 
like magnetic anisotropy and giant magnetoresistence (GMR) are much larger in low dimen- 
sional systems. Whereas the magnetic anisotropy is in general larger in two dimensional 
systems due to the reduced symmetry and thus a different crystal potential, the change of 
the magnetic moments is generated by the magneto- volume effect, i.e. the dependence of 
the magnetic moment on the atomic volume |]I^|-|TH]. The latter is affected by the differ- 



ences between the equilibrium lattice constants of substrate and overlayer and the induced 
overlayer structures. One of the most striking examples is the Fe/Cu(001) system which 
shows a rich magnetic phase diagram in the range from 1 to 11 monolayers |ll7|-[T9|l. Therein 
heavily distorted fee structures appear. The interlayer distances are relaxed (resulting in 
a fct structure) and shifts of the atomic positions in and perpendicular to the layer plane 
occur pOf. Also the in-plane lattice constants of consecutive layers are relaxed. For the 
Fe/Cu(001) system experiments confirmed the sensitivity of NOLIMOKE on structural 
changes of the top layer. 

Since all these effects take place in configurations with only a few monolayers it is pos- 
sible to obtain them directly by ab initio methods. A lot of calculations were done for 
freestanding and supported monolayers, most of them aiming at the computation of mag- 



netic anisotropy. For 3d transition metals the effects of hybridization with the substrate [21 



different coordination and rf-band filling [^,^, the stability of the monolayer for different 



magnetic configurations, such as ferromagnetic, antiferromagnetic or canted spin |24-26 



and superstructures |23] were calculated from first principles. 

Since the theoretical prediction of the sensitivity of nonlinear optics to surface mag- 
netism []28|-pO[] and the first experiments the applicability of NOLIMOKE to thin film 
systems has been demonstrated for several phenomena. Due to the fact that SHG is also 
generated at buried interfaces, properties of different interfaces in multilayer systems could 
be separated pT|J5^ . For that purpose it was important that the size of the nonlinear 
magneto-optical effects, namely the nonlinear Kerr rotation, is strongly increased compared 
to linear optics These measurements also show a dependence on interface rough- 

ness By direct comparison of linear and nonlinear MOKE, changes of the magnetic 

properties of the topmost layer during the growth process were detected for the Co/Cu(001) 
system [^, since NOLIMOKE is sensitive to the surface and interlayer only while linear 
MOKE integrates over the magnetism of all layers. Also quantum well states which occur in 
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sandwich structures could be resolved. This has been shown both experimentally [0,p%-|39 



and theoretically [0], by exploiting the fact that NOLIMOKE spectra reflect characteris- 
tic features of the bandstructures Recent work predicted that even the influence of 



surface antiferromagnetism on the optical signal can be resolved by the nonlinear magneto- 
optical Kerr-effect (NOLIMOKE) |^. This has already been shown before for SHG of the 
antiferromagnetic non-inversion symmetric bulk Cr203 . 

One further potential of SHG, which to our knowledge has not been applied to magnetic 
systems so far, is the strongly enhanced sensitivity to submonolayer coverages [^,|^. Second 
harmonic generation by small particles is enhanced by local-field effects. In the case of 



clusters deposited on a substrate this gives rise to signals for particle sizes around 1 nm pS 



which is far beyond the resolution limit of linear optics. For spherical particles the effects of 
local-field enhancement are well known by the linear Mie theory . Extensions to nonlinear 
optics show an enhanced sensitivity of the size- dependent resonances compared to the 
linear case In the case of 3d transition-metal overlayers it should be possible to resolve 
nanostructures of nm size with low density by making use of the submonolayer coverage 
sensitivity of SHG and the different in-plane symmetries of the nanostructures and the 
substrate. From the experimental point of view the preparation of nanostructures can now 
be achieved by state of the art techniques such as molecular beam epitaxy varying the growth 
parameters (e.g. the deposition rate or the temperature [|9,50|). 

So far calculations of SHG generated by metal surfaces are mainly restricted to simple 
and noble metals which are well described by the model of a free electron gas. These systems 
were intensely studied by Liebsch and coworkers [pl| -[53[|. They also calculated anisotropic 
contributions [^^ -|56| and the influence of steps ||57|-p9|] and obtain good agreement with 
experiment. Other authors studied the change of the SHG yield in the presence of adsorbates 
on simple metal surfaces within density functional theory pD|-p^. For these nonmagnetic 
systems the intraband transitions show stronger contributions than the interband transition. 
Thus a better model for the screening effects is necessary, whereas in the case of transition 
metals the response is mainly due to interband transitions. Then the intraband effects can 
be added by applying a Drude model using experimental parameters . Calculations of the 
linear magneto-optical Kerr-effect (MOKE) indicate that ab initio methods including spin- 
orbit coupling and an highly accurate determination of the dipole transition matrix elements 
are necessary to obtain magneto-optical spectra which can be compared to experimental 
values [|6^ -|65[|. To some extent this was realized for nonlinear magneto-optics by Pustogowa 
et al. In their work the Kerr spectra of Fe films with one to seven layers and the 

dependence of the Kerr-spectra of a Fe(OOl) monolayer on the in-plane lattice constant have 
been calculated by determining the electronic bands within a full potential linear muffin tin 
orbital (FP-LMTO) code. Spin-orbit coupling was treated within first-order perturbation 
theory and the optical matrix elements were approximated as constants. 

Here we will use the FLAPW method and go beyond this work within two respects: differ- 
ent orientations of Fe-monolayers ((001), (110), (111) of fee) are investigated and, apart from 
the Kerr-spectra and the magnetic moments, we focus on the spin-orbit coupling strength 
and its structural dependence. Though microscopically both spin-orbit coupling and spin- 
polarization are necessary to generate magneto-optical response, spin-orbit coupling plays a 
special role, since the spin-orbit coupling strength is directly proportional to the size of the 
magneto-optical Kerr effect. This is known from studies of linear MOKE Thus e.g. 
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the Kerr-rotation of an Fe/Pt system is much larger than that of an Fe layer, since the large 
spin-orbit coupling of Pt contributes via hybridization with the magnetic Fe-layer |]64,69 



This knowledge is important for applications in storage technology, where magneto-optics 
is applied in a configuration, where a perpendicular easy axis in combination with an in- 
creased Kerr-rotation is preferred. In contrast to the spin-orbit coupling the dependence on 
the magnetic moment is rather complicated. Nevertheless little is known about the spin-orbit 
coupling constants of thin film systems contrary to their magnetic moments. 

In our work the optical spectra are determined by using the same approximations as 
in i.e. the matrix elements are taken as constants and the effects of spin-orbit coupling 
in the wavefunctions are treated within first-order perturbation theory. Since the spin-orbit 
induced changes of the wavefunctions yield first-order effects |^ while spin-orbit induced 
shifts of the eigenenergies give rise to second-order effects, we neglect spin-orbit coupling 
in the calculations of the electronic bands, which are obtained from first principles. The 
validity of this approach will be shown below. 

By the choice of the investigated monolayers, we want to study several aspects of struc- 
tural changes. Firstly, we investigate the influence of relaxation of the in-plane lattice 
constant, which is varied over a wide range for the Fe(OOl) monolayer. Secondly, substrates 
of different orientations are simulated by comparing the results for the Fe(OOl), Fe(llO) 
and Fe(lll) monolayers, which also reveals the role of coordination. These structures are 
deduced from the bulk fee lattice. Two lattice constants are considered, the lattice constant 
induced by Cu fee bulk and an even smaller value. Thirdly, the role of nanostructuring is 
studied for regular arrays of stripes, which can be created by viewing the closed monolayer 
gular array of chains and then relaxing the distance between the chains. Although 
this structure is rather artificial, it reveals the effect of reducing the dimension of the layer 
in a second direction. Also we compare our results with previous calculations. 

Future work will address the calculation of the optical dipole matrix elements to get 
the full information on the size of the NOLIMOKE spectra and to exploit the symmetry 
properties of the systems, which will be of special interest in the calculation of special nanos- 
tructures like triangular islands. This includes the determination of the lateral resolution 
limit of nonlinear optics. 

The paper is organized as follows: In Sec. II, we will outline the theory for the nonlinear 
magneto-optical response and our method to calculate the spin-orbit coupling. Then the 
result part follows, which is divided into three subsections (Sees. Ill A, B, C), each for 
the comparison of different characteristic changes of the structures. The paper ends with 
summary and outlook (Sec. IV). 
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II. THEORY 



Within the electric-dipole approximation the polarization P of the medium can be ex- 
panded in terms of the incident field E as 

= x'i^E, + x^lE.Ek + ... , 

where x^^^ ^ind x*^^"* are the linear and second harmonic susceptibilities We calculate the 
nonlinear magneto-optical response within the theoretical framework introduced by Hiibner 
and Bennemann p9| and obtain the nonlinear susceptibility in the electric dipole approxi- 
mation as 



— ie^ 



(2q, 200) = -— {(k + 2q, /" |z|k, /) (k, /|j|k + q, /') (k + q, l'\k\k + 2q, /") 

k,/,/',«" 

/(gk+2q.i")-/(^k+q.i') _ f(Ey,+ ^,l')-fiE^,l) 

^k+2q,/" - ^ki - 2huj + i2hai 



The indices z, j, k run over x, y and z. In previous calculations the wave-functions and 
the band energies were calculated neglecting spin-orbit coupling. Instead spin-orbit coupling 
was taken as a perturbation and the product of the three matrix elements were calculated 
using first order perturbation theory to yield 

^ (k + 2q, I) (k, /|j|k + q, /') (k + q, r|A;|k + 2q, /") (2) 

where the wavefunctions and energies do not contain spin-orbit coupling and the spin-orbit 
coupling constant is taken from the atomic value of the spin-polarized ci-bands. The matrix 
elements are approximated as constants. This approach includes explicit inversion symmetry 
breaking but makes it impossible to distinguish the different elements of the tensor Xijk- 
Nevertheless the resulting nonlinear susceptibility 

/(-^k+2q|| ,l"<T)-/(-^k+q|| ,1'ct) /(-^k+q|| ,!'<t )^ /(-E^k,i<T ) 

^(2)^2q 2UJ M) - H ^k+2^ll,l"--^k+'Hl,l'''-^+'^°'^ ~ -gk+q||.i'^--gk.ia-?^^+^ 

^ cr k,i,F,F' -Ek+2q|| ,/"cr — -EkZcr — 2/iu; -I- z2/;,ai 

(3) 

reflects the spectral dependence of a magnetic tensor element, since spin-orbit coupling 
enters in first order. Nonmagnetic tensor elements (and all even order tensor elements) also 
consist of the zeroth order (and the corresponding higher even orders) in spin-orbit coupling. 
Thus they do not contribute to magneto-optics within first order and yield larger values. 
Due to our approximations we add in Eq. (^ a spin index a, drop the indices which specify 
the tensor elements and add the factor originated by the approximate size of the matrix 
elements. 

The susceptibility is exclusively built on interband transitions. We will use this approx- 
imation throughout this paper, since interband resonances dominate the optical response of 
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metallic systems. Thus in our case we will call the dependence of LjHmx^'^\uj) on the photon 
energy NOLIMOKE "spectra" . For details we refer to . 

Calculations including spin-orbit coupling will only affect the band-energies E\^, because 
the factor Aso describes the effect of SOC in the wavefunctions and the matrix elements, 
which are included in x*^^) as constants C, are not calculated explicitly. 

In this work the bandstructures are obtained from first principles using the full potential 
linear augmented plane wave (FLAPW) method WIEN95 |jl2|. Additionally we implemented 
spin-orbit coupling in a second variational step as described e.g. by |pTi|7Tl , |72| . After the 
self-consistent determination of the wavefunctions and eigenenergies (quantities which are 
obtained self-consistently without SOC are marked by a suffix "sc" in the following) the 
Hamiltonian matrix is determined including spin-orbit coupling 



ISC I i sc 



ij ij 

to obtain the eigenfunctions 

= EPn(?)'^k„ g=l,2, ... 

n 

and the corresponding eigenenergies e(g) shifted by spin-orbit coupling (g is the index of the 
eigenenergies and Pn{(l) is the coefficient of the n-th basis-function in the g-th eigenfunction) . 
Here, spin-orbit coupling is not calculated self consistently, especially the basis functions are 
not affected by SOC. The procedure is known to yield good agreement with exact results W3 



To determine the spin-orbit part of the Hamiltonian the basis functions of the FLAPW 
method have to be taken into account. The basis set consists of the standard basis functions 



E [Aim (k,) ui (r. El) + Bim (k^) in (r, Ei)] Yi^ r < R^t 

Im (4) 



and the so called local orbitals, which are introduced to describe the low lying semi-core 
states irgj?! 



E [Aim (ki) ui (r. El) + Bim (kO in (r, Ei) + Q„ (k^) in (r, E2)] Yim r < Rmt 

Zm . (5) 

r > Rmt 



and are included for all /-values for which semi-core states appear (-Rmt is the muffin tin 
radius). The radial functions are obtained from the Schrodinger equation 



2 a _ ^ ((i + i) ^ 



<(r) = Elulir) , (6) 



where the localization energies E" are chosen to be at the center of the band. 
The spin-orbit operator 



Hso = ■ (VV X p 



2 

[a is the fine structure constant) is applied in the spherical approximation: 
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r or 

dV 

— = 0, r > RrrA , 

or 

since the gradient of the potential yield its largest contributions near the core, where the 
potential is almost spherical. This yields 

, ,ldV ^IdV 
Hso = — s • (r X p) -— = — s ■ L- — . 

2 r or 2 r or 

The spin-orbit matrix elements 



«2 ^ 1 dV\ 



are calculated by separating the angular and radial parts. This yields 

= E { ALA*m(ki)Am'(k,-) + 

[AL (k.) Bi^, (k,) + fi*„ (k,) A,^, (k,)]} 
I YC^ (f) s • l.Yim' (r) 



Imm' 



with the spin-orbit coupling constants 



\l = 

uu — 


2 


/ 

"'r<Rmt 


dr 
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(r) 




\l = 

uu — 
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dr 




(r) 
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Kil = 
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dr 


«r 


(r) 


^ 5^ (^) 
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\/ = 

uu — 
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^uu = 
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/ 



(8) 



and the spin-orbit coupling constants including local orbital functions: 

dr uf (r) r^^uf (r) 

dV I I 
dr uf (r) r—^^uf (r) 

dV I I 
dr ill (r) r—^^iif (r) 

Thus we get three (six when local orbitals are involved) spin-orbit coupling constants for one 
Z- value which are formed by a radial integral over the radial part of the basis functions and 



(9) 
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the radial derivative of the potentiaL Furthermore one has to take into account that spin- 
orbit couphng mixes the spins, thus the Hamilton-matrix gets off-diagonal elements within 
the space of the spin | and | basis functions. Spin is not a good quantum number anymore 
and the wavefunctions consist of both spin | and spin | contributions. The spin-orbit matrix 
elements Eq. get then additional spin indices o and a' . 




Therein the spin-index of the potential is equal to the spin index of the basis function on 
the right, since the spin-orbit operator acts on it. The fact that the potentials are different 
for the spins, but the basis functions are not, leads to the requirement to make the matrix 
explicitly hermitian, since the spin-orbit operator is. This affects only the spin-flip matrix 
elements. 
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III. RESULTS 



In Sec. III. A we will simulate the effect of lattice relaxation. This can be achieved 
experimentally by different substrates, assuming pseudomorphic growth. We will show NO- 
LIMOKE spectra of free-standing Fe(OOl) monolayers with in-plane lattice constants varied 
from a = 2.4 A, which is slightly below the value of the nearest-neighbor distance in Cu 
fee bulk a=2.56 A, to a = 2.76 A, which is close to the nearest-neighbor distance of Fe 
bcc bulk. For the comparison of the trends of the magnetic moments and the spin-orbit 
coupling constants we extend the range of lattice constants from 2.22 A to 3.18 A, the latter 
corresponding approximately to the value of bcc W. In Sec. III.B the same quantities are 
shown for Fe monolayers with different structures, i.e. the fee (111), (001) and (110), which 
are schematically displayed in Fig. |I[ The structures are studied for the Cu fee nearest- 
neighbor distance a=2.56 A and a = 2.4 A . It should be possible to get a measure of the 
structural changes from the NOLIMOKE spectra. In Sec. III.C we will show the influence 
of nanostructuring on the NOLIMOKE spectra by analyzing stripe structures as indicated 
in Fig. H. The Fe(lll) monolayer can be interpreted as an array of "zig-zag" -stripes. To 
reduce the dimension of the structure we vary the distance of the stripes, which is indicated 
by d, where in the case oi d = h the layer is equal to the (111) structure. 



A. Fe(OOl) monolayers 

Fig. 1^ shows the NOLIMOKE spectra cu^Imx*^^^ [uj] of the Fe(OOl) monolayer as a function 
on the in-plane lattice constant. As analyzed by Pustogowa et al. within a tight binding 
scheme the first maximum and the zero are mainly due to features of the d bands, whereas 
for higher photon energies the role of the s-p-bands is more dominant. In particular they 
showed that the position of the zero is a measure for the ci-band width and the height of the 
maximum is proportional to the magnetic moment. Thus the d bands generate the features 
of the spectra in the optical region. In our case the zeros show a clear dependence on the 
lattice constant. The positions shift to lower energies with increasing lattice constant. Since 
this point characterizes the (i-bandwidth, the bandwidth is reduced upon lattice expansion. 
From the bandstructure it can also be seen that bands above the range of visible frequencies 
are shifted to lower energies with increasing lattice constant, which generates the different 
slopes in the high-energy part of the spectra. The height of the maxima starts to increase 
with the spin-polarization for lattice constants from 2.4 A to 2.58 A . For the larger lattice 
constants (a=2.67 and 2.76 A) there are no more significant changes of the peak height, 
as can be seen in the inset of Fig. ^ However, the position of the peaks is shifted in 
proportion to the value of the magnetic moments to lower energy values. The dependence of 



the maximum on the magnetic moments agrees with previous works |67,68|, where no clear 



dependence of the linear magneto-optical response on the size of the magnetic moments was 



found. Additionally, from a tight binding calculation, Pustogowa et al. |^5[ found a linear 
dependence of the maximum for magnetic moments between and 2.5 /i^, but a similar 
behavior for moments between 2.5 /is and 3.4 ^b- The difference should reflect that in both 
calculations the magnetic moments are changed by different mechanisms. Whereas in the 
tight binding calculations the magnetic moments were affected by changes of the exchange 
coupling constant J, in our case the magnetic moments are varied by changing the lattice 
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constant, which does not only shift the relative positions of the d subbands, but also their 
width. 

The values of the magnetic moments increase with increasing lattice constant. This is 
shown in Fig. ^ where the size of the magnetic moments (filled circles) is plotted as a 
function of the in-plane lattice constant in units of fis- If interpolated our results agree 



very well with calculations by Wang et al. [Q, who obtained 3.04yUB for a lattice constant 



of a = 2.56 A, and with results by Gay and Richter [^], who obtained 3.20/iB for a lattice 
constant of a = 2.88 A compared to our results of 3.08 for a = 2.56 A and 3.24/iB for 
a = 2.88 A, respectively. Nevertheless one has to keep in mind that the error in the magnetic 
moments is around ±5% due to the chosen accuracy in our calculations. 

In Fig. ^ we also compare the values of the magnetic moments directly with the spin- 
orbit coupling constants A^"^^ t| and || defined in Eq. (|]), i.e. the spin-orbit coupling 
constants for bands with d character for | and | spin. These constants are the important 
ones for magneto-optics since the (i-bands exhibit the magnetic moments. Combinations 
of the radial functions other than {ui,ui) are of less interest since the radial dependence 
of wavefunctions is mainly described by the ^/-functions. The plot shows two important 
properties of the coupling constants: (i) The values increase with decreasing lattice constants 
and (ii) the difference between the coupling constants of | and | spin show a clear dependence 
on the magnetic moments. Inspection of the potentials for the different lattice constants 
shows that the size of the spin-orbit coupling constants is not directly governed by changes 
of the potential, i.e. the derivative of the potential shows no changes near the core, where 
the largest values of the derivative occur. As a consequence, the changes of the constants 
must be induced by changes of the radial functions. This is shown in Fig. |^, where the 
square of the function M;(r) and the integrand uidV/drui defined in Eq. (P) are plotted as 
a function of the radial distance. The insets show that the increase of the maximum of 
the integrand, which causes the changes of the coupling constants with decreasing lattice 
constants, is proportional to the changes of the square of ui (the maxima of both quantities 
are normalized to one). Thus the changes of the coupling constants are caused by changes 
of the potential near the muffin tin radius Rmt? which alters the probability of the maximum 
of the radial functions also close to the nuclei. In addition the dependence of the differences 
between the |t and J,| coupling constants on the magnetic moments reveals the variation of 
the potential with different spin-subband occupation and by changing the radial functions 
via Eq. (|]). For even larger values of the lattice constants the coupling constants should 
reach the atomic value, which is approximately 50 meV. In the case of Ge, the spin-orbit 
splitting of the Ap electrons in the solid is 0.43 eV at the F point, a 30% enhancement 
to the spin-orbit splitting of 0.21 eV in the Ge atom. By comparing the wavefunctions and 
potentials in the solid and the atom, we find that this increase in the spin-orbit coupling 
strength in the solid in Ge is caused by a quite different reason. Because of the covalent 
bond the charge distribution is not only enhanced between the atoms but also near the 



core. Daalderop et al obtained the coupling constants of bulk Fe and find a much 
larger difference between the coupling constants for | and J, spin and also the values differ 
slightly. Since they used a LMTO code the differences should mainly be due to the different 
definitions of the coupling constants resulting from the different basis sets used. 

In Tab. | the values of the additional spin-orbit coupling constants for / = 1 and 2 and 
within the combinations of the radial functions {ui, ui) and (m;, iii) are listed. The dependence 
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of the coupling constants with / = 1 on the lattice constant differs significantly from the 
values for / = 2. Their changes are much more pronounced, namely the values decrease by 
about 40% rather than by only 5% for / = 2. Also a spin-polarization appears only for the 
largest value of the lattice constants whereas it slightly changes for I = 2. These properties 
reveal that the p states are much more influenced by the changed binding characteristics. 
Though in the case of Fe monolayers the p band is not occupied and its center is tens 
of eV above the Fermi-level, p-states could become accessible by optical excitations via 
hybridization with low lying s- and p-bands of an appropriate substrate such as Mo or W. 
In these systems the large values of the SOC-constants could, in particular in the case of 
pump-probe femtosecond spin dynamic experiments, significantly affect spin-orbit induced 
spin-flip contributions by strong excitations. 

In principle contributions related to the radial functions iii should not contribute sig- 
nificantly since the values of the coefficients Bim are in general much smaller than Aim, 
nevertheless the changes of the coupling constants A^^ should reflect some features of the 
shape of the bands. Since the values for / = 2 show no spin-polarization, the shape of the 
subbands should be nearly equal, also the derivative increases quite strongly with the lattice 
constants indicating narrower bands. For / = 1 the increase is even stronger in agreement 
with the values for A^"^^, the smaller values compared to Z = 2 reflect stronger dispersed 
bands. 

A direct comparison of the spectra of the Fe(OOl) monolayer for a = 2.76 A obtained 
within the FLMTO method and our FLAPW method in Fig. |^ shows good agreement in 
the region of low photon energies. The position of the maximum is near 1.5 eV in both cases 
and the energy where the susceptibility crosses zero is 3 eV. In both calculations the same 
model for the nonlinear magneto-optical susceptibility was applied. Thus the differences in 
the region of higher photon energies should be an effect of the different ab initio methods 
and in this special case due to the different basis sets. In the LMTO-method the number of 
basis functions is much smaller than in the FLAPW method which leads to a lack of bands 
high above the Fermi-level unless the calculations are performed for several localization 
energies. This is in agreement with the fact that for spectra which are based on the Fe 
bulk bandstructure we find no significant differences in both methods also for high photon 
energies. 

Our results obtained for the changed lattice constants agree very well with the results by 
Pustogowa et al. for the same system. Since they calculated spectra for lattice constants 
larger than 2.76 A their changes of the zero are smaller due to the nonlinear dependence of 
the shape of the bands on the lattice constants, which is also refiected by the dependence of 
the magnetic moments on the lattice constants in Fig. ^. In contrast to their calculations in 
our case the position of the maximum shows a clear dependence on the lattice constant. 

The optical spectra also depend on the type of approximations applied to the calculations 
of the electronic bands. This can be seen in Fig. 0, where the spectrum of a Fe(OOl) 
monolayer with a = 2.76 A is calculated using different approximations for the exchange 
correlation potential. We compare the generalized gradient approximation (GGA), which is 



used for all calculations throughout this work, in the parameterization by Perdew et al. ||80 



with the local spin density approximation (LSDA) in the parameterization of Perdew and 



Wang 1 81]. Since GGA corrects for overbinding, the bandwidth should be lowered and thus 



the zero should be at lower energy. In our case the opposite behavior occurs. The LSDA 
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values are lower in energy. The spectra obtained with the different LSDA approximations 
show no strong deviations. GGA yields a slightly higher magnetic moment which can be 
responsible for the higher value of the maximum and the larger value of the zero. In general 
it is expected that spin-orbit coupling counteracts GGA, since the bandwidth increases by 
spin-orbit induced shifts. But in the case of the NOLIMOKE spectra the effects of spin-orbit 
coupling which enters the spectra via the bandstructure are negligible as can be seen in Fig. 
where the NOLIMOKE spectrum of Fe(llO) is plotted both with and without SOC effects 
on the bandstructure. This reflects that spin-orbit induced changes of the eigenenergies only 
contribute in second order to the spectrum [ ]7D| . 



B. Fe(OOl), (110) und (111) monolayers 



In Figs. I and |TO|the NOLIMOKE spectra of the Fe(OOl), (110) and (111) monolayers 
are compared for the nearest-neighbor distance of Cu fee bulk, a=2.56 A, and for a=2.4 A . 
Since the nearest-neighbor distance is equal in the different structures, the changes reveal 
the different coordination, which is six in the hexagonal (111) layer, four in the square lattice 
(001) and two in the rectangular lattice (110). The different coordination determines the 
area of the two dimensional (see Fig. |I|) unit cell containing one atom to ^f2c? for the (110) 
structure, c? for (001) and \/2)/2a^ for (111). The next nearest neighbor distance is \/2a in 
the (110), 2a in the (001) and VSa in the (111) layer. 

For both nearest-neighbor distances it can be seen that the lattice with the lowest coordi- 
nation shows the smoothest spectra, whereas for highest coordination the most complicated 
structure appears. This is a general aspect of coordination and can also be seen in the 
bandstructures [^. Between eV and the zero around 3 to 4 eV the (110) spectrum has 



a sinusoidal shape, in the (001) spectra first a maximum followed by a shoulder appears, 
whereas in the (111) case a dominant maximum is surrounded by two shoulders. The dif- 
ferences are much more pronounced for the smaller nearest-neighbor distance, in the case of 
a=2.56 A the spectra are closer, as can be seen e.g. at the zero point, and it is more difficult 
to define a maximum. Comparing the positions of the zero for both lattice constants one 
notices that the shifts are larger for higher coordination. Whereas the zero remains more 
or less constant in the case of Fe(llO), it is shifted to lower energies by approx. 0.2 eV for 
Fe(OOl) and 0.4 eV for Fe(lll). Thus one can say that the dependence of the spectra on 
the lattice constants is proportional to the coordination. 

Roughly the same holds for the magnetic moments in Figs. |TT| and There the magnetic 
moments and the spin-orbit coupling constants ]] and [[ are plotted for three 
different coordination numbers corresponding to (111), (001) and (110) layers. The magnetic 
moments of the Fe(llO) layers change by only approx. 0.1 hb compared to approx. 0.25 hb 
for the (111) and (001) monolayers. As expected the values of the magnetic moments 
increase with lowered coordination. Comparison of the spin-orbit coupling constants show 
that the changes induced by the different coordination are quite small compared to the 
changes induced by different nearest-neighbor distances. Thus one can say that in a first 
approximation the values of the spin-orbit coupling constants depend on the nearest-neighbor 
distance and remain constant for different coordination. The values of the coupling constants 
with / = 1 and the combination of the radial functions (u;, iii) confirm this since they show 
no significant changes with the coordination as can be seen from Tabs. |T| and |T|. 
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C. Stripe structures 



The NOLIMOKE spectra of the stripe structures of Fig. ^ are plotted in Fig. [T^ for 
different distances d of the stripes in comparison to the spectrum of the Fe(lll) monolayer, 
which can be interpreted as a stripe structure with distance d = h (see Fig. H). The spectra 
show no behavior which can be easily interpreted in terms of the bandwidth, corresponding 
to a zero point in the spectra, or the magnetic moment, corresponding to a maximum at 
a certain position in the low-energy-regime. Only the more complicated structure of the 
spectra for d > h is clear from the lifting of degeneracies in the bandstructure, which results 
from the breaking of symmetry. Thus, the number of bands increases for d > h, since 
there are two nonequivalent atoms in the unit cell, compared to one for Fe(lll). Since the 
bandstructures of the different stripe structures show no pronounced overall changes, the 
differences in the spectra should be an effect of the details of the bands. 

If one first neglects the spectra for the stripes with the largest distance (long-dashed line) 
the behavior is quite regular in the sense that the maximum value of the spectra decreases 
with increasing distance and that the zero point shifts to lower energies and reaches zero for 
d = h + 1.04 A . By drastically increasing the "interstripe" distance to d = h + 1.78 A the 
NOLIMOKE spectra differ from these trends and exhibit a shape which is similar to the 
spectra of the closed layers in Fig. 0, and |TU|. This can be interpreted as an oscillatory 



behavior of the electronic structure with the distance. The change of the trends for the 
distance oi d = h + 1.78 A can also be observed in Fig. |T3| for the spin-orbit coupling 
constants and more or less also for the magnetic moments. Since the difference of their 
values for d = h + 1.04 A and = -f 1.78 A is very small, fi remains essentially constant. 

The relatively small changes of the spin-orbit coupling constants with I = 2 shown in 
Fig. |l^ and Tab. |^ imply the same interpretation as in Sec. III.B for the layers with dif- 
ferent coordination. In a first approximation the values of the coupling constants remain 
unchanged and thus their values depend mainly on the nearest-neighbor distance, which is 
fixed here by the constant structure of the isolated stripes. Thus, for tailoring the SOC- 
constants and the magnetic moments, the choice of the substrate will be much more efficient 
than nanostructuring while keeping the nearest-neighbor distance constant. Nevertheless 
nanostructuring can still have a strong effect on dynamical properties of the spin and mag- 
netic moments, for which the size of the SOC-constants is quite important. For certain 
island sizes a phase transition from ferro- to superparagnetism will occur and affect the spin 
dynamics. 

The behavior of the coupling constants with / = 1, which are also listed in Tab. |IV|, show 
one difference. The values are also nearly constant for distances larger than h, but they 
show a clear increase, when the distance changes from d = h to d = h + 0.46 A, i.e. at the 
onset of the reduction of the symmetry. Thus the SOC-constants of the p-bands are not 
only much larger but also much more sensitive to the bond characteristics than the (i-bands. 

Clearly further investigations are necessary to understand the strong changes of the 
NOLIMOKE spectra for the different stripe distances. 
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IV. SUMMARY AND OUTLOOK 



We presented results for the structural dependence of the nonlinear magneto-optics of Fe 
monolayers. Since in our theory the optical dipole-transition matrix elements are approxi- 
mated as constants, we cannot analyze the symmetry properties of the nonlinear suscepti- 
bility tensor. Instead we focused on the spectral dependence of a magnetic tensor clement, 
the magnetic moments, and the spin-orbit coupling constants, the latter two reflecting the 
microscopic origin of magneto-optics. 

In the case of the Fe(OOl) monolayer-spectrum the characteristic features like the position 
of the first zero, which is related to the d-band width, and the position of the maxima are 
shifted to lower energies with decreasing lattice constant. The changes are stronger for 
smaller lattice constants, which also holds for the magnetic moments and the spin-orbit 
coupling constants. The values of the maxima, which should be related to the magnetic 
moments show no clear trend. The spectra of layers with different coordination numbers 
show characteristic differences in the shape as well as in the position of the maximum and 
the first zero. The differences are more pronounced for smaller lattice constants. Reducing 
the dimensionality of the monolayers, simulated by onedimensional stripes with different 
"interstripe" -distances, results in dramatic changes of the spectra. Their shapes show no 
similarities with the monolayer spectra any more. In contrast the values of the spin-orbit 
coupling constants depend in a first approximation only on the nearest-neighbor distance. 
This was shown for both the Fe monolayers with different coordination number and the 
onedimensional stripe structures. 

For the Fe(OOl) monolayers the SOC-constants show the opposite behavior as the mag- 
netic moments, they increase with decreasing lattice constants. As an important result the 
difference between the coupling constants for | and J, spin is proportional to the magnetic 
moments, caused by the dependence of the potentials on the occupation of the subbands. 
The values of the magnetic moments show results well-known for itinerant ferromagnets. 
Increasing the lattice constants or decreasing the coordination enhances their values. The 
same holds for increasing the distances of stripes in the quasi-onedimensional structures. 

Our results clearly show the strong dependence of the NOLIMOKE spectra on structural 
changes and also indicate that the spectral dependence of the magneto-optical response is a 
valuable source of information on the structure of the investigated system. 

Future work will address the completion of our ab initio theory by the optical transition 
dipole matrix elements, which is of major importance not only for the determination of ab- 
solute signal values but also for the study of structural dependencies, since the susceptibility 
tensor reflects the symmetry of the system. Also we will investigate structures with further 
reduced dimensionality, i.e. zerodimensional islands, and apply nonlinear magneto-optics to 
antiferromagncts. 

We acknowledge financial support by Deutsche Forschungsgemeinschaft through Sfb 290 
and by TMR Network NOMOKE contract No. FMRX-CT96-0015. 
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FIGURES 

FIG. 1. Geometry of the different layers investigated in this paper. The twodimensional unit 
cell containing one atom is shown, deduced from the fee lattice in the (110), (001) and (111) 
direction. 

FIG. 2. Chains built by stretching the Fe(lll) monolayer as indicated. The unit cell used 
for the bandstructure calculation is indicated by the solid rectangle. To simulate the chains the 
distance d is increased compared to h. In the case of d = h the layer is equal to the Fe(lll) 
monolayer. 

FIG. 3. NOLIMOKE spectra of Fe(OOl) monolayers with the lattice constant varying from 

a = 2.4 A to a = 2.76 A . 

FIG. 4. Magnetic moments and spin-orbit coupling constants for Fe(OOl) monolayers as a 

function of the lattice constant. up and X''^ dn denote the SOC-constants obtained from the 

matrixelements within the spin-combinations || and ||. 

FIG. 5. Square of the radial function U2{r) and integrand U2{r)'^ • dV/dr as a function of the 

radius. In both cases the highest values are normalized to unity. The scale of the y-axes in the 

insets are equal, which indicates that changes in the integrand are directly caused by the radial 

functions. 

FIG. 6. NOLIMOKE spectra in the case of a Fe(OOl) monolayer obtained by Pustogowa et 
al. [75] using an FP-LMTO code and by the present authors using the FLAPW method. 

FIG. 7. Comparison of the NOLIMOKE spectra of a Fe(OOl) monolayer, obtained by using 
GGA [76], LSDA with the parameterization of [77]. 

FIG. 8. Comparison of the NOLIMOKE spectra obtained for the Fe(llO) monolayer without 

spin-orbit coupling and with band shifts induced by SOC. The effect of SOC in the wavefunctions 

via the optical dipole matrix elements has not been taken into consideration. 

FIG. 9. NOLIMOKE spectra as a function of the photon energy of the fundamental light for 

Fe(OOl), Fe(llO) and Fe(lll) monolayers using the Cu fee lattice constant a=2.56 A . 

FiG. 10. NOLIMOKE spectra as a function of the photon energy of the fundamental light for 

Fe(OOl), Fe(llO) and Fe(lll) monolayers using the lattice constant a=2.4 A . 

FIG. 11. Magnetic moments and spin-orbit coupling constants X^^ for the spin-combinations 
tt and II of the Fe(llO), Fe(OOl) and Fe(lll) monolayers with the nearest-neighbor distance of 
Cu fee bulk a = 2.56 A . The values are plotted as a function of the coordination. 

FIG. 12. Magnetic moments and spin-orbit coupling constants for the spin-combinations 

II and II of the Fe(llO), Fe(OOl) and Fe(lll) monolayers with the nearest-neighbor distance 

a = 2.4 A . The values are plotted as a function of the coordination. 

FIG. 13. NOLIMOKE spectrum of the different stripe structures for distances between the 

stripes varying from /i to /i -|- 1.78 A as described in Fig. ^. 

FIG. 14. Magnetic moment and spin-orbit coupling constant X^^ for the spin-combinations || 
and II for the different stripe structures as a function of their "interstripe" distance d. 
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TABLES 



TABLE I. Values of the spin-orbit coupling constants AJ^„ and A^^ with Z = 1,2 for || and J, J, 
spin-combinations for the Fe(OOl) monolayer as a function of the lattice constants a between 2.4 
and 2.76 A . 

TABLE IL Values of the spin-orbit coupling constants AJ^„ and A^^ with Z = 1,2 for || and 
11 spin-combinations for the Fe(lll), Fe(OOl) and Fe(llO) monolayer with the Cu fee bulk near- 
est-neighbor distance a = 2.56 A . 

TABLE in. Values of the spin-orbit coupling constants A^^ and A^^ with I = 1,2 for || and 
II spin-combinations for the Fe(lll), Fe(OOl) and Fe(llO) monolayer with the nearest-neighbor 
distance a = 2.4 A . 

TABLE IV. Values for the spin-orbit coupling constants AJj„ and A'^,^ with I = 1,2 for || and 
II spin-combinations of the stripe structures as a function of the "interstripe" distance d. 
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Tab. 1 



[eV] 


a=2.40 A 


a=2.49 A 


a=2.58 A 


a=2.67 A 


a=2.76 A 




\(— 1 


TT 


0.40789 


0.33736 


0.29471 


0.26022 


0.29714 




\l=l 


ii 


0.41127 


0.33583 


0.29551 


0.26304 


0.31008 




\l=l 

iiu 


TT 


0.00082 


0.00018 


0.00061 


0.00128 


0.00601 




\l=l 

uii 


ii 


0.00115 


0.00041 


0.00102 


0.00186 


0.00697 




\l=2 
Ku 


TT 


0.06182 


0.06039 


0.05941 


0.05865 


0.05909 




\l=2 
Ku 


ii 


0.06321 


0.06193 


0.06140 


0.06091 


0.06095 




\l=2 


TT 


0.01372 


0.01604 


0.01868 


0.02164 


0.02459 






ii 


0.01385 


0.01621 


0.01883 


0.02175 


0.02489 





Tab. 2 



[eV] 


Fe(lll) 


Fc(OOl) 


Fc(llO) 




^1=1 II 


0.34339 


0.34822 


0.34601 




\/=l 1 1 


0.35045 


0.35664 


0.35519 




ALf TT 


0.00244 


0.00262 


0.00252 




\i=i 1 1 


0.00299 


0.00322 


0.00316 




^uu TT 


0.05989 


0.06058 


0.06023 




\l=2 1 1 

uu 4- 4- 


0.06191 


0.06220 


0.06207 




^uii TT 


0.01807 


0.01784 


0.01790 




\l=2 1 1 


0.01818 


0.01802 


0.01802 




Tab. 3 








[eV] 


Fe(lll) 


Fe(OOl) 


Fe(llO) 




^^uu TT 


0.40746 


0.40789 


0.40900 




\/=i 1 1 

"^uu 4- 4- 


0.41047 


0.41127 


0.41360 




^uiL TT 


0.00084 


0.00082 


0.00086 




\/=l 1 1 

uii 4-4- 


0.00114 


0.00115 


0.00123 




TT 


0.06175 


0.06183 


0.06176 




\l=2 1 1 

uu 4- 4- 


0.06305 


0.06321 


0.06350 




Auit TT 


0.01380 


0.01372 


0.01370 




\l=2 1 1 
uii 4- 4- 


0.01392 


0.01385 


0.01378 
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Tab. 4 



[eV] 




d^h + 0.46 A 


d^h + 1.04 A 


d=/i + 1.78 A 






0.35653 


0.38331 


0.38857 


0.38502 




\l=l 1 1 


0.35388 


0.38784 


0.39421 


0.39004 




AL"ii^ TT 


0.00008 


0.00127 


0.00146 


0.00130 




\i=i 1 1 


0.00022 


0.00166 


0.00189 


0.00171 




Ar.' TT 


0.06063 


0.06105 


0.06120 


0.06085 




UU 


0.06261 


0.06277 


0.06290 


0.06239 




ALf TT 


0.01518 


0.01507 


0.01502 


0.01504 




\l=2 1 1 


0.01523 


0.01519 


0.01514 


0.01519q 
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